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Abstract 

Background: Herbivore coexistence on the Late Cretaceous island continent of Laramidia has been a topic of great 
interest, stemming from the paradoxically high diversity and biomass of these animals in relation to the relatively 
small landmass available to them. Various hypotheses have been advanced to account for these facts, of which 
niche partitioning is among the most frequently invoked. However, despite its wide acceptance, this hypothesis has 
not been rigorously tested. This study uses the fossil assemblage from the Dinosaur Park Formation of Alberta as a 
model to investigate whether niche partitioning facilitated herbivorous dinosaur coexistence on Laramidia. 
Specifically, the question of feeding height stratification is examined in light of the role it plays in facilitating 
modern ungulate coexistence. 

Results: Most herbivorous dinosaur species from the Dinosaur Park Formation were restricted to feeding no higher 
than approximately 1 m above the ground. There is minimal evidence for feeding height partitioning at this level, 
with ceratopsids capable of feeding slightly higher than ankylosaurs, but the ecological significance of this is 
ambiguous. Hadrosaurids were uniquely capable of feeding up to 2 m quadrupedally, or up to 5 m bipedally. There 
is no evidence for either feeding height stratification within any of these clades, or for change in these ecological 
relationships through the approximately 1.5 Ma record of the Dinosaur Park Formation. 

Conclusions: Although we cannot reject the possibility, we find no good evidence that feeding height 
stratification, as revealed by reconstructed maximum feeding heights, played an important role in facilitating niche 
partitioning among the herbivorous dinosaurs of Laramidia. Most browsing pressure was concentrated in the herb 
layer, although hadrosaurids were capable of reaching shrubs and low-growing trees that were out of reach from 
ceratopsids, ankylosaurs, and other small herbivores, effectively dividing the herbivores in terms of relative 
abundance. Sympatric hadrosaurids may have avoided competing with one another by feeding differentially using 
bipedal and quadrupedal postures. These ecological relationships evidently proved to be evolutionarily stable 
because they characterize the herbivore assemblage of the Dinosaur Park Formation through time. If niche 
partitioning served to facilitate the rich diversity of these animals, it may have been achieved by other means in 
addition to feeding height stratification. Consideration of other feeding height proxies, including dental microwear 
and skull morphology, may help to alleviate problems of underdetermination identified here. 
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Background 

During the Late Cretaceous, a shallow inland sea divided 
North America into two longitudinally- arrayed land- 
masses. The eastern landmass— called Appalachia [1] — 
supported an enigmatic fauna known only from scant 
fossil remains [2]. The western landmass— called 
Laramidia [1]— supported a rich diversity of herbivorous 
dinosaurs. Included among these were such small forms 
(< 100 kg) as hypsilophodontids, pachycephalosaurids, and 
leptoceratopsids. Putatively facultative herbivores (or omni- 
vores) included troodontids, oviraptorosaurs, ornithomimids, 
and therizinosaurids [3]. The megaherbivorous dinosaurs— 
primarily represented by ankylosaurs, ceratopsids, and 
hadrosaurids— form the most significant component of this 
fauna [4-6], in terms of both population (75-82% of total 
dinosaur fauna; [7]) and body size. 

This high diversity of contemporaneous herbivores in 
Laramidia has long puzzled investigators. Laramidia had 
a total estimated area between just 4 million km 2 [8] and 
7.7 million km 2 [5], and dinosaur distribution was segre- 
gated into distinct northern and southern faunal prov- 
inces that presumably reflect strict habitat preferences 
linked to a palaeoclimatic gradient [4-6,9]. By contrast, 
the area of sub-Saharan Africa, where living vertebrate 
herbivores are most diverse, is approximately 23.6 
million km 2 , with the highest diversity occurring in the 
East African savannah [10]. The diverse herbivorous 
dinosaur fauna of Laramidia therefore appears to have 
been spatially restricted compared to modern analogs, 
and this likely would have increased the potential for 
resource competition. The problem of the coexistence of 
these dinosaurs is further exaggerated when considering 
the presumably large nutritional requirements of the 
megaherbivorous forms [11-14], as well as their high 
population densities [15-22], both of which would have 
placed increased pressure on the resource base. 

Two main hypotheses have been proposed to account 
for the coexistence of so many herbivores in Laramidia. 
The first hypothesis contends that plant resources were 
simply not limiting, a scenario that could be explained by 
numerous causal factors. For example, it may have been 
that, in spite of the fact that the nutritional requirements 
of the megaherbivorous dinosaurs were absolutely high, 
the inferred bradymetabolic thermoregulatory systems of 
these animals imparted relatively low nutritional require- 
ments compared to their mammalian counterparts, 
minimizing pressure on the resource base [5,14,23]. Alter- 
natively, plant resources may not have been limiting be- 
cause primary productivity was elevated during the Late 
Cretaceous [23,24]. This suggestion was also advanced to 
account for the increased abundance of browsers during 
the Miocene [25,26]. Increased primary productivity may 
have resulted in part from the prevalence of marginal 
coastal environments during the Late Cretaceous [27]; 



however, palaeoclimatological modeling [28,29] and 
experimental evidence [30] also suggest that elevated at- 
mospheric pressure and C0 2 during the Late Cretaceous 
would have yielded similar effects. A third option might 
be that predation pressure from tyrannosaurids and 
dromaeosaurids was sufficiently high during the Late 
Cretaceous to suppress herbivore population densities, 
leading to reduced pressure on the resource base as in the 
first scenario. A similar mechanism is thought to help 
shape the structure of modern African ungulate commu- 
nities [31]. 

The second hypothesis maintains that plant resources 
were limiting, and that herbivorous dinosaur coexistence 
was facilitated by dietary niche partitioning [6,13,32]. If so, 
these animals should have been adapted in such a way as to 
avoid competing with one another, perhaps varying in 
morphology so as to specialize on different plant types 
[7,8,33-42]. Numerous ecomorphological correlates 
among the herbivores have been hypothesized to reflect 
niche partitioning, including differences in snout shape 
[7,33,34,36-41,43,44], jaw mechanics and tooth shape 
[36,40,45,46], and feeding height stratification [12,13,33,34]. 

The late Campanian-aged Dinosaur Park Formation 
(DPF) of Alberta preserves the most diverse dinosaur 
assemblage currently known, including the most diverse 
assemblage of large-bodied herbivores from Laramidia 
[47]. Currently, approximately 33 herbivorous and om- 
nivorous dinosaur species are recognized [48]. Biostrati- 
graphic work shows that these taxa are distributed 
heterogeneously throughout the formation, such that 
different species are restricted to different horizons 
within the formation [47,49-51]. Although not all 33 
taxa co-existed, over half likely co-existed in certain time 
intervals, including at least eight megaherbivorous taxa 
(Megaherbivore Assemblage Zone 1 [51]). Therefore, 
this assemblage provides a unique opportunity to 
evaluate possible dietary niche partitioning among Late 
Cretaceous herbivorous dinosaurs. 

Niche partitioning via feeding height stratification is 
thought to have facilitated sauropod diversity during the 
Late Jurassic [33,34,44,52-58]. Likewise, some [12,13,33] 
have suggested that the rich diversity of herbivorous di- 
nosaurs in the Late Cretaceous of Western Interior 
North America was fostered, in part, by a similar mech- 
anism, but rigorous quantification of these patterns and 
formal tests of specific hypotheses are lacking. The aim 
of this paper is to test the hypothesis that the long-term 
coexistence of these animals was facilitated by dietary 
niche partitioning, more specifically via feeding height 
stratification. 

Methods 

We examined the question of feeding height stratifica- 
tion by estimating the maximum feeding heights (MFHs) 
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of the dinosaurs from the DPF. It is prudent to note 
from the outset that MFH does not necessarily reflect 
habitual feeding height. Habitual feeding height is very 
difficult to determine in fossil forms, particularly in light 
of the fact that amniotes rarely maintain the cervical col- 
umn in an osteologically neutral pose [59]. Regardless, 
most discussions to date regarding the feeding habits of 
dinosaurs from the DPF have focused on MFH, which is 
our primary concern. See Discussion below for further 
treatment of maximal and habitual feeding heights. 

The data used in this study are listed in Additional file 1. 
We excluded juvenile specimens, identified by their small 
size and/or under-developed cranial ornamentation, be- 
cause we were specifically interested in estimating the 
MFHs of fully-grown individuals. Presumably, juveniles fed 
at intermediate heights and their inclusion in this study 
would have no effect on the calculation of MFHs. We also 
excluded the facultatively herbivorous/omnivorous thero- 
pods because they formed only a small component of the 
DPF fauna [7,12,60,61], and were therefore unlikely to have 
consumed an appreciable amount of plant matter. Many of 
these animals are also known only from fragmentary 
remains, and their MFHs can be estimated only with a 
great deal of uncertainty. 

We studied all fossil specimens with permission from 
the following institutions: ACM, Beneski Museum of 
Natural History, Amherst College, Amherst; AMNH, 
American Museum of Natural History, New York; 
CMN, Canadian Museum of Nature, Ottawa; MCSNM, 
Museo Civico di Storia Naturale di Milano, Milan; NSM 
PV, National Museum of Nature and Science, Tokyo; 
ROM, Royal Ontario Museum, Toronto; TMP, Royal 
Tyrrell Museum of Palaeontology, Drumheller; UALVP, 
University of Alberta Laboratory of Vertebrate Palaeon- 
tology, Edmonton; YPM, Yale Peabody Museum, New 
Haven. 

Maximum feeding height estimation 

When estimating the MFHs of the herbivorous dino- 
saurs, it was important to devise some metric that 
allowed for the use of both articulated and disarticulated 
skeletons, thereby maximizing sample size. MFH can 
often be measured directly from mounted specimens, 
but many reasonably complete skeletons remain 
disarticulated in museum collections. For this reason, we 
calculated MFH as follows: for the quadrupedal forms, 
we estimated MFH from shoulder height, calculated by 
adding the lengths of the humerus, radius, and metacar- 
pal III. Observations of mounted skeletons reveal that 
this is an appropriate proxy for MFH because the mouth 
is approximately level with the glenoid of the shoulder 
(Figure 1A, B). In the case of ceratopsids, the mouth 
could not be brought much higher than this because the 
extension of the head at the atlanto-occipital joint would 




Figure 1 Quadrupedal MFH proxies used in this study. Shoulder 
height (calculated as the combined lengths of the humerus, radius, 
and metacarpal III) is used to approximate the MFHs of ankylosaurs 
(A), ceratopsids (B), and other small, quadrupedal forms. Hip height 
(calculated as the combined lengths of the femur, tibia/fibula, and 
metatarsal III) is used to approximate the quadrupedal MFHs of 
hadrosaurids (C). Dashed lines indicate the approximate quadrupedal 
MFH. Image (B) provided by G. Danis (used with permission). 



eventually cause the huge parietosquamosal frill to abut 
the dorsal surface of the shoulders (Figure IB). For 
hadrosaurids, quadrupedal MFH was estimated from hip 
height, calculated by summing the lengths of the femur, 
tibia, and metatarsal III (the fibula was used if the tibia 
was not available). Again, observations of mounted 
hadrosaurid skeletons confirm that this is an appropriate 
proxy for quadrupedal MFH (Figure 1C). 

For bipedal forms (including the facultatively bipedal 
hadrosaurids and leptoceratopsid), bipedal MFH was cal- 
culated using trigonometry: assuming that the animal in 
question could rear up on its hind legs only so far as the 
distal tip of its tail would touch the ground (a reasonable 
assumption, given that the tails of many of these animals 
were relatively inflexible due the presence of either 
tendon trellises or ossified myorhabdoi [62] that ex- 
tended much of their lengths), a right triangle is formed 
(Figure 2) with the vertical side equal to hip height and 
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Figure 2 Depiction of the trigonometric model used here to 
estimate bipedal MFHs of hadrosaurids. See Methods for details. 



the hypotenuse equal to tail length. The intersection of the 
hypotenuse and the horizontal side of the triangle forms 
angle 6, which can be calculated as 0 = sin" 1 (hip height/tail 
length). The triangle can then be extended by adding the 
combined length of the trunk and neck to the hypotenuse 
to give the full body length. The product of sin 6 and body 
length yields an estimate for bipedal MFH. 

We took measurements from the literature when they 
could not be made directly. Examples of complete, artic- 
ulated vertebral columns from the DPF are nonexistent. 
For this reason, we assumed that mounts that included 
restored vertebrae were reasonably accurate and were 
included in this study; however, entirely missing or 
restored sections of the vertebral column (e.g., cervicals, 
thoracics, caudals) or limb elements were not used. We 
estimated these missing elements from more complete 
specimens using reduced major axis regression (RMA; 
Additional file 2). RMA is preferable to ordinary least 
squares (OLS) regression because it does not assume that 
the independent variable is measured without error. There 
is some debate, however, about whether RMA is preferable 
to OLS when extrapolating beyond the dataset [63], as we 
occasionally did here. For species lacking postcranial ma- 
terial, we modeled MFHs after closely related taxa of simi- 
lar size (e.g., Unescoceratops koppelhusae was modeled 
after Leptoceratops gracilis [CMN 8889], "Stegoceras" 
("Prenocephale") breve and the unnamed pachycephalosaur 
after S. validum [UALVP 2], and cf. Orodromeus after 
Parksosaurus warreni [ROM 802]). 

The various sources of error associated with the 
reconstruction of MFHs must be noted. For instance, 
habitual limb flexion in any of the taxa considered here 
would result in overestimated MFHs. Ceratopsids, in 
particular, are generally thought to have held their fore- 
limbs in a semi-sprawled posture [64-68], which would 



result in smaller MFHs. Conversely, our inability to 
account for the influence of epiphyseal cartilage on feed- 
ing height calculations would result in underestimated 
MFHs. Holliday et al. [69] recently showed that, among 
archosaurs, epiphyseal cartilage significantly increases 
limb length relative to that calculated from the bones 
alone, the results differing by up to 10%. By accounting 
for missing cartilage in multi-tonned dinosaurs using vari- 
ous "cartilage correction factors" (CCFs), Holliday et al. 
[69] found that the heights of these animals may be con- 
sistently underestimated by up to 0.5 m. Unfortunately, it 
is not possible to determine exactly how much epiphyseal 
cartilage is missing from every individual, and to apply dif- 
ferent CCFs to separate taxa in the reconstruction of max- 
imum feeding heights would tell us more about the 
influence of CCFs in producing (or reducing) statistical 
differences, and less about how the fossil bones them- 
selves differ. On the other hand, it is unlikely that the 
same CCF should apply across all taxa. In the best case 
scenario, the confounding effects of limb flexion and 
epiphyseal cartilage on the calculation of maximum 
feeding heights would cancel each other out. Overall, 
the error in reconstructing MFH is likely on the 
decimetre scale. Although this may affect the results 
in some cases (e.g., discerning between low-level 
browsers), we believe that this level of error should 
permit the major patterns to be discerned. 

Statistical comparisons 

We compared MFHs at coarse (family/suborder), medium 
(subfamily/family), and fine (genus) taxonomic scales. We 
did not consider the species level because sample size was 
generally too low at this resolution to permit meaningful 
statistical comparisons. We used the non-parametric 
Kruskal-Wallis test because samples were typically quite 
small (n < 30) and non-normal. This test lacks the power of 
parametric tests, and is more prone to committing Type II 
errors (reporting false negatives), but is more robust to 
committing Type I errors (reporting false positives). We 
conducted posthoc pairwise comparisons using Mann- 
Whitney U tests with Bonferroni correction. Bonferroni 
correction was designed to counteract the problem of mul- 
tiple comparisons, wherein the probability of committing 
a Type I error increases with the number of simultan- 
eous comparisons being made [70]. This problem is 
rectified by multiplying the j^-value by the number of 
pairwise comparisons, effectively lowering the signifi- 
cance level. However, because Bonferroni correction 
provides little power and is probably too conservative 
[70,71], we also report uncorrected probabilities for the 
purpose of interpretation. We performed all statistical 
procedures using the software program PAST 2.12 [72]. 

Because the DPF does not represent a single assem- 
blage of contemporaneous organisms, time-averaging is 
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an issue. This has the effect of masking palaeoecological 
patterns that are otherwise distinguishable only at fine 
temporal resolutions [73]. For this reason, we minimized 
the effects of time-averaging by making the above 
comparisons within each of the two most inclusive 
Megaherbivore Assemblage Zones (MAZs) identified by 
Mallon et al. [51]. To summarize, MAZ-1 encompasses 
the lower 28 m of the DPF, and MAZ-2 encompasses 
intervals from 29-52 m. Although this time-constrained 
approach theoretically increases the probability of recov- 
ering differences that would otherwise be masked by the 
effects of time-averaging, there is a trade-off in that 
sample size (and hence statistical power) is reduced con- 
siderably. Also, this approach does not completely re- 
move the effects of time-averaging because the MAZs 
are themselves time-averaged over a period of approxi- 
mately 600 Ka [51]. Although we took steps to minimize 
the amount of time-averaging in this study, the remaining 
time bias is likely too large to effectively capture true, 
interacting palaeocommunities. It is possible to reduce the 
time bias further by dividing the MAZs into approxi- 
mately 300 Ka sub-zones [51], more closely approxi- 
mating true palaeocommunities, but the sample size 
(and resulting statistical power) per sub-zone becomes 
drastically reduced in doing so. 



Results 

Time-averaged approach 

The MFHs for the time-averaged analysis are depicted in 
Figure 3A. The most striking aspect of this diagram is 
that, irrespective of body mass, most species were 
restricted to feeding below approximately 1 m from the 
ground. The leptoceratopsid Unescoceratops appears to 
have been particularly limited in its feeding range be- 
cause it was unable to reach above 0.5 m while in a 
quadrupedal posture. However, this animal probably 
could have reared bipedally to reach somewhat beyond 
1 m in a manner similar to other coexisting, small or- 
nithischians (cf. Orodromeus and pachycephalosaurids). 

With respect to the megaherbivorous forms, the aver- 
age quadrupedal MFHs of ankylosaurs generally fall 
slightly below 1 m, whereas those of ceratopsids plot 
slightly above 1 m. Hadrosaurids clearly have the most 
distinct MFHs, which consistently plot above 2 m in a 
quadrupedal posture, or 4 m in a bipedal posture. 
Hadrosaurines (Gryposaurus and Prosaurolophus) have 
slightly taller average bipedal MFHs than lambeosaurines, 
but the uncertainty intervals for all hadrosaurids in this 
posture are quite large and overlap substantially. This may 
reflect variation in reconstructed portions of the spinal 
column and estimated missing elements. 

The Kruskal-Wallis test reveals highly significant 
differences among all higher-level taxa with adequate 
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Figure 3 Reconstructed MFHs of the herbivorous dinosaurs 
from the DPF. A, time-averaged analysis; B, MAZ-1 analysis; 
C, MAZ-2 analysis. Scale bars represent 95% confidence intervals. 



representation (N = 62, H = 55.82, p < 0.001); Ankylosauria, 
Ceratopsidae, and Hadrosauridae each differ significantly 
from one another in quadrupedal MFH (Table 1). As 
expected, the differences are exacerbated when bipedal 
hadrosaurids are considered, which differ significantly from 
all quadrupedal postures. No further significant differences 
are recovered with increasingly finer taxonomic reso- 
lution, even between the apparently different bipedal 
hadrosaurines and lambeosaurines (Tables 2-3). 
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Table 1 Coarse-scale Mann-Whitney pairwise comparisons of the time-averaged analysis 



Ankylosauria Ceratopsidae Hadrosauridae (Q) Hadrosauridae (B) 



Ankylosauria 




0.00 


0.00 


0.00 


Ceratopsidae 


0.00 




0.00 


0.00 


Hadrosauridae (Q) 


0.00 


0.00 




0.00 


Hadrosauridae (B) 


0.00 


0.00 


0.00 





N = 62, H = 55.82, p = 4.554 x 1 0" 12 . 

Bonferroni corrected p-values shown in lower left triangle; uncorrected pairwise comparison p-values shown in upper right triangle. Significant results reported in 
bold. Abbreviations: Q, quadrupedal posture; B, bipedal posture. 



Time-constrained approach 

The results of the time-constrained analyses (MAZ-1 
and -2) largely mirror those of the time-averaged analysis; 
the same MFH relationships are maintained (Figure 3B, C). 
Note that some taxa, while known from the DPF, are not 
represented in either MAZ-1 or -2 due either to a lack of 
appropriate material or imprecise provenance data (e.g., 
ankylosaurs in MAZ-2, Chasmosaurus in MAZ-1). Prob- 
abilities are reduced in the time-constrained analyses, often 
below significant levels, due to correspondingly reduced 
sample sizes (Tables 4, 5, 6, 7, 8 and 9). 

Discussion 

Feeding height stratification in extant ecosystems 

Ecological separation of coexisting species is achieved 
along the axes of food, time, and space [74,75]. For 
example, the African savannah biome supports 31 species 
of large (> 5 kg), herbivorous mammals [76], and competi- 
tion among its members is typically alleviated via the 
selection of different food types, the occupation of the 
same habitat at different times, or the occupation of differ- 
ent habitats at the same time [77-79]. The ecological 
separation of these ungulates may also be achieved along 
a vertical gradient, with different species feeding at differ- 
ent heights within the canopy or grass cover [77-83]. For 
example, Bell [78] demonstrated that, among grazers, 



Table 2 Medium-scale Mann-Whitney pairwise 
comparisons of the time-averaged analysis 





An 


No 


Ce 


Ch 


Ha (Q) 


Ha (B) 


La (Q) 


La (B) 


An 




0.51 


0.00 


0.01 


0.00 


0.00 


0.00 


0.00 


No 


1 




0.11 


0.19 


0.01 


0.01 


0.00 


0.00 


Ce 


0.04 


1 




0.67 


0.00 


0.00 


0.00 


0.00 


Ch 


0.30 


1 


1 




0.01 


0.01 


0.00 


0.00 


Ha (Q) 


0.10 


0.40 


0.07 


0.40 




0.00 


0.36 


0.00 


Ha (B) 


0.06 


0.30 


0.04 


0.30 


0.10 




0.00 


0.30 


La (Q) 


0.01 


0.11 


0.01 


0.11 


1 


0.01 




0.00 


La (B) 


0.01 


0.11 


0.01 


0.11 


0.02 


1 


0.00 




N = 62, H 


= 56.17, p = 8, 


.721 X10 


10 











Bonferroni corrected p-values shown in lower left triangle; uncorrected 
pairwise comparison p-values shown in upper right triangle. Significant results 
reported in bold. Taxonomic abbreviations: An, Ankylosauridae; Ce, 
Centrosaurinae; Ch, Chasmosaurinae; Ha, Hadrosaurinae; La, Lambeosaurinae; 
No, Nodosauridae. Other abbreviations as in Table 1. 



zebra (Equus burchelli) tend to ingest the tallest, most fi- 
brous portions of the herb layer, wildebeest {Connochaetes 
taurinus) and topi (Damaliscus korrigum) select the more 
nutritious middle layer, and Thomsons gazelle (Gazella 
thomsoni) take in fruits from the ground. With the excep- 
tion of the topi, this grazing succession is reflected by 
the decreasing body size of the animals (but see [83]). 
Similarly, feeding height stratification has been said to 
operate among browsing ungulates [77,79,81]. This hypoth- 
esis was tested explicitly on a subset of African browsers by 
du Toit [82], who found that giraffe (Giraffa cameloparda- 
lis), kudu (Tragelaphus strepsiceras), impala (Aepyceros 
melampus), and steenbok (Raphicerus campesteris) are 
stratified by mean feeding-height. On average, giraffe feed 
at heights between 2 and 3 m, kudu at heights near 1 m, 
and impala and steenbok at heights below 0.5 m. Although 
there was significant overlap between the browsing heights 
of some of these species, du Toit [82] observed that feeding 
height stratification is more pronounced during the dry sea- 
son, when the use of woody browse is increased to com- 
pensate for the reduced availability of green forage in the 
herb layer. The long neck of the giraffe has been cited as an 
example of an adaptation to escape competition occurring 
at lower browsing heights [84]. 

Feeding height stratification in the Dinosaur Park 
Formation assemblage 

As documented in extant taxa, numerous studies have 
commented on the importance of feeding height stratifica- 
tion as a mechanism for niche partitioning among herbiv- 
orous dinosaurs. For example, within the Upper Jurassic 
Morrison Formation of the western United States, up to 
five sauropod genera are thought to have lived in sympatry 
[85]. Feeding height stratification has been repeatedly in- 
voked as a means of facilitating their coexistence, with dif- 
ferent sauropods using their long necks to feed at 
different heights within the environment. Evidence for this 
interpretation includes reconstructed neck morphology 
[33,34,56,57], tooth wear analysis [44,53-55,58,86], and 
jaw mechanics [52,54]. 

Feeding height stratification has likewise been invoked 
to account for the diverse herbivorous dinosaur fauna of 
the DPF [12,13,33], but no formal test of this hypothesis 
has been conducted to date. Notably, most treatments of 
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Table 3 Fine-scale Mann-Whitney pairwise comparisons of the time-averaged analysis 





Euo 


Pan Cen Sty Cha Gry (Q) 


Gry (B) 


Pro (Q) 


Pro (B) 


Cor (Q) 


Cor (B) 


Lam (Q) 


Lam (B) 


Par (Q) 


Par (B) 


Euo 




0.59 0.01 0.07 0.03 0.01 


0.01 


0.07 


0.03 


0.00 


0.00 


0.01 


0.01 


0.07 


0.07 


Pan 


1 


0.39 0.11 0.38 0.03 


0.03 


0.11 


0.05 


0.01 


0.01 


0.03 


0.03 


0.11 


0.11 


Cen 


0.85 


1 0.33 0.55 0.02 


0.02 


0.08 


0.04 


0.01 


0.01 


0.02 


0.02 


0.08 


0.08 


Sty 


1 


1 1 0.77 0.11 


0.11 


0.25 


0.15 


0.06 


0.06 


0.11 


0.11 


0.25 


0.25 


Cha 


1 


1 1 1 0.05 


0.05 


0.15 


0.08 


0.02 


0.02 


0.05 


0.05 


0.15 


0.15 


Gry (Q) 


1 


1111 


0.03 


0.82 


0.05 


0.92 


0.01 


0.67 


0.03 


0.82 


0.11 


Gry (B) 


1 


1111 1 




0.11 


0.60 


0.01 


0.51 


0.03 


0.89 


0.11 


0.49 


Pro (Q) 


1 


1111 1 


1 




0.15 


0.19 


0.06 


0.11 


0.11 


0.70 


0.25 


Pro (B) 


1 


1111 1 




1 




0.02 


0.36 


0.05 


0.86 


0.15 


0.77 


Cor (Q) 


0.36 


1 0.61 1 1 1 




1 






0.00 


0.51 


0.01 


0.31 


0.06 


Cor (B) 


0.36 


1 0.61 1 1 1 




1 




0.23 




0.01 


0.92 


0.06 


0.38 


Lam (Q) 


1 


1111 1 




1 




1 


1 




0.03 


0.25 


0.11 


Lam (B) 


1 


1111 1 




1 




1 


1 


1 




0.11 


0.82 


Par (Q) 


1 


1111 1 




1 




1 


1 


1 


1 




0.25 


Par (B) 


1 


1111 1 




1 




1 


1 


1 


1 


1 




N = 59, H = 53.52, p = 1.553 x10~ 6 . 

Bonferroni corrected p-values shown in lower left triangle; uncorrected pairwise comparison p-values shown in upper right triangle. Significant results reported in 



bold. Taxonomic abbreviations: Cen, Centrosaurus; Cha, Chasmosaurus; Cor, Corythosaurus; Euo, Euoplocephalus; Gry, Gryposaurus; Lam, Lambeosaurus; Pan, 
Panoplosaurus; Par, Parasaurolophus; Pro, Prosaurolophus; Sty, Styracosaurus. Other abbreviations as in Table 1. 



feeding height stratification make little or no mention of 
the small ornithischians that inhabited the Late Cretaceous 
landscape of Laramidia. This is probably because these 
forms are both poorly known and relatively rare in the fos- 
sil record. However, Brown et al [48] recently demon- 
strated a substantial taphonomic size bias in the record of 
the DPF, and it is therefore likely that small ornithischians 
comprised a larger proportion of the herbivore fauna 
than previously assumed. The MFHs of these animals 
were probably restricted to less than approximately 1 m, 
which may have placed them in competition with the 
megaherbivores discussed below. 

Ankylosaurs 

Beland and Russell [12] and Coe et al [13] suggested that 
ankylosaurs from the DPF exhibited different MFHs, with 
Euoplocephalus feeding on herbs below 0.5 m and 
Panoplosaurus feeding on woody vegetation up to 1 m. 
However, the current data do not support this hypothesis. 
The mean MFH of Euoplocephalus is closer to 0.8 m, 



whereas that of Panoplosaurus approximates 0.9 m. There 
is also considerable overlap of MFH between these taxa, 
with some specimens of Euoplocephalus reaching as high 
as 0.88 m, and some specimens of Panoplosaurus only 
reaching 0.77 m. Apart from the reduced ankylosaur sam- 
ple (see below), this moderate overlap of MFHs likely ac- 
counts for the fact that Euoplocephalus and Panoplosaurus 
do not differ significantly from one another. The conten- 
tion of Weishampel and Norman [87] that ankylosaurs 
were generally restricted to browsing between 1 and 2 m is 
not supported here. By extending the forelimb proportions 
of Euoplocephalus (AMNH 5403) to the 542 mm-long 
humerus of Ankylosaurus (AMNH 5214), the largest 
known ankylosaur, the total forelimb length is estimated to 
be only 1.1 m. 

Ceratopsids 

Several authors [12,13,87,88] have likewise suggested 
that ceratopsids were capable of reaching heights of 2 m; 
however, no evidence was provided for this value. 



Table 4 Coarse-scale Mann-Whitney pairwise comparisons of the MAZ-1 analysis 





Ankylosauria 


Ceratopsidae 


Hadrosauridae (Q) 


Hadrosauridae (B) 


Ankylosauria 




0.24 


0.00 


0.00 


Ceratopsidae 


1 




0.00 


0.00 


Hadrosauridae (Q) 


0.00 


0.02 




0.00 


Hadrosauridae (B) 


0.00 


0.02 


0.00 





N = 36, H = 31.1 1,p = 8.049 x 10~ 7 . 

Bonferroni corrected p-values shown in lower left triangle; uncorrected pairwise comparison p-values shown in upper right triangle. Significant results reported in 
bold. Abbreviations as in Table 1. 
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Table 5 Medium-scale Mann-Whitney pairwise 
comparisons of the MAZ-1 analysis 





An 


No 


Ce 


Ha (Q) 


Ha (B) 


La (Q) 


La (B) 


An 




0.82 


0.11 


0.15 


0.15 


0.04 


0.04 


No 


1 




0.67 


0.05 


0.05 


0.01 


0.01 


Ce 


1 


1 




0.05 


0.05 


0.01 


0.01 


Ha (Q) 


1 


1 


1 




0.08 


0.55 


0.01 


Ha (B) 


1 


1 


1 


1 




0.01 


0.80 


La (Q) 


0.87 


0.12 


0.12 


1 


0.30 




0.00 


La (B) 


0.87 


0.12 


0.12 


0.30 


1 


0.00 




N = 36, H = 31 .25, p = 2.269 x 1 0" 5 . 

Bonferroni corrected p-values shown in lower left triangle 


uncorrected 





pairwise comparison p-values shown in upper right triangle. Significant results 
reported in bold. Abbreviations as in Tables 1 and 2. 

Instead, our results indicate that ceratopsids more likely 
browsed no higher than approximately 1 m, as suggested 
by Dodson [7]. One well-preserved specimen of Tricera- 
tops (NSM PV 20379), among the largest ceratopsids that 
ever existed, could not reach above 1.2 m, either [68]. 

Hadrosaurids 

The common claim that hadrosaurids could reach heights 
up to 4 m [7,12,13,33,87,88] is supported by our results. 
The largest hadrosaurids from the DPF, Gryposaurus and 
Prosaurolophus, probably could reach heights approaching 
5 m in a bipedal posture (Figure 3). Despite these maxima, 
some authors [87,89] have proposed that hadrosaurid 
feeding was probably concentrated below 2 m, which 
would accord with the quadrupedal feeding postures 
calculated here. 

Both vertebrate microfossil and skeletal remains sug- 
gest that hadrosaurids formed approximately 40% of the 
herbivorous dinosaur assemblage in the DPF ecosystem 
[7,61], which likely translates to a greater proportion of 
the herbivore biomass because hadrosaurids are the 
largest members of the fauna. The remainder of herbivores, 



notably the large ceratopsids and ankylosaurs which com- 
bined form approximately an equivalent proportion in 
terms of relative abundance, fed at or below the 1 m mark. 
Although feeding heights could not be discriminated below 
1 m, the fact that hadrosaurids could reach up to 5 m and 
were therefore segregated from all other herbivores is likely 
fundamental in partitioning the resource base. Importantly, 
hadrosaurids were capable of reaching shrubs and low- 
growing trees that were beyond the reach of ceratopsids, 
ankylosaurs, and other small herbivores, effectively dividing 
the herbivores in terms of relative abundance. This may 
also have allowed hadrosaurids to escape resource stresses 
imparted by low browsers, and may have facilitated the co- 
existence of large herds of ceratopsids and highly abundant 
hadrosaurids [20,90] in DPF palaeoecosystems. 

Dinosaur browsing and vegetation structure 

The prevailing climate of the DPF has been described as 
warm temperate, as revealed by tree-, tooth-, and bone- 
growth ring data [91,92], sedimentological data [93], and 
biogeographic data [94,95]. Regional leaf physiognomic 
data have also reinforced this interpretation [96]. This 
climate is, in part, thought to have given rise to both 
open and closed habitats [96-99] akin to those of mod- 
ern ecosystems [36]. Braman and Koppelhus ([100]:124) 
describe the landscape of the DPF as having been "wet 
everywhere, at least for portions of the year", with dense 
vegetation lining the rivers, and more open habitats 
occurring further distally. 

The open habitats of the DPF and surrounding regions 
were likely dominated by ferns and low-growing angio- 
sperms. Coe et al. ([13]: 235) even proposed the exist- 
ence of extensive "fern prairies", analogous to modern 
grasslands, but Tiffney [101] stressed that evidence for 
such fern-dominated communities is lacking. Nonethe- 
less, Wing et al. [102] subsequently reported one excep- 
tional fossil flora from the mid-Maastrichtian in which 



Table 6 Fines-scale Mann-Whitney pairwise comparisons of the MAZ-1 analysis 





Pan Cen Gry (Q) 


Gry (B) 


Cor (Q) 


Cor (B) 


Lam (Q) 


Lam (B) 


Par (Q) 


Par (B) 


Pan 


0.67 0.05 


0.05 


0.01 


0.01 


0.11 


0.11 


0.11 


0.11 


Cen 


1 0.05 


0.05 


0.01 


0.01 


0.11 


0.11 


0.11 


0.11 


Gry (Q) 


1 1 


0.08 


0.90 


0.03 


0.39 


0.15 


0.77 


0.15 


Gry (B) 


1 1 1 




0.03 


0.90 


0.15 


0.77 


0.15 


0.77 


Cor (Q) 


0.64 0.64 1 






0.01 


0.62 


0.07 


0.40 


0.07 


Cor (B) 


0.64 0.64 1 




0.23 




0.07 


0.87 


0.07 


0.40 


Lam (Q) 


1 1 1 




1 


1 




0.25 


0.70 


0.25 


Lam (B) 


1 1 1 




1 


1 


1 




0.25 


0.70 


Par (Q) 


1 1 1 




1 


1 


1 


1 




0.25 


Par (B) 


1 1 1 




1 


1 


1 


1 


1 




N = 34, H = 29.26, p = 5.86 x 1 0" 4 . 

Bonferroni corrected p-values shown in lower left triangl 


^; uncorrected pairwise comparison p-values shown in upper right triangle. Significant results reported in 



bold. Abbrevations as in Tables 1 and 3. 
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Table 7 Coarse-scale Mann-Whitney pairwise comparisons of the MAZ-2 analysis 

Ceratopsidae Hadrosauridae (Q) Hadrosauridae (B) 

Ceratopsidae 0.03 0.01 

Hadrosauridae (Q) 0.08 0.05 

Hadrosauridae (B) 0.04 0.16 

N = 13, H = 10.38, p = 5.559 x10~ 3 . 

Bonferroni corrected p-values shown in lower left triangle; uncorrected pairwise comparison p-values shown in upper right triangle. Significant results reported in 
bold. Abbreviations as in Table 1. 



ferns and other "pteridophytes" account for nearly 50% 
of the total ground cover. By comparison, these same 
plants account for approximately 40% of the total 
palynomorph abundance in the DPF [103]. 

The Late Cretaceous saw the radiation of the angio- 
sperms, which typically took the form of "weedy" herbs 
and shrubs growing in open or marginal habitats 
[12,33,97,98,101,104,105]. Angiosperm trees, although 
inferred to have existed elsewhere [96], probably did not 
occur in the DPF, as evidenced by the lack of diagnostic 
fossil wood [100]. It is commonly argued ([98]: 125 
[104,106]) that angiosperms occurred most regularly in 
coastal and fluvial depositional settings, occupying 
"stream-side and aquatic habitats, the forest understory and 
early successional thickets"; however, Wheeler and Lehman 
[107] noted the existence of angiosperm-dominated com- 
munities in southern upland environments as well, where 
conifers were otherwise thought to have dominated 
[98,104]. By virtue of their r-selected life history strategies, 
it is likely that angiosperms were capable of growing in a 
wide variety of habitats [101]. 

A common theme in the literature is the persistence of 
open-habitat cycadophytes (bennettitaleans and cycads) 
as forage for Late Cretaceous herbivores [34,97,108-113]. 
However, bennettitaleans went extinct by the Santonian 
[114], and cycads were probably absent from the DPF 
([65,115,116], D. R. Braman, pers. comm., 2012), having 
been replaced by angiosperms [98,117]. It is by no 
means clear that cycadophytes were common enough 
elsewhere during the Late Cretaceous to support dense 
megaherbivore populations, but their seed coats may 
have served to supplement dinosaur diets [97,109-111]. 



Unlike tropical forests, temperate forests typically 
exhibit limited stratification [118], and there is little rea- 
son to suspect that the temperate forests of the DPF 
were any different. Wolfe and Upchurch [96] proposed 
that such forests were, in fact, relatively sparse, with 
sunlight often penetrating fully through to the ground. 
Palynofloral and macroplant evidence from the DPF sug- 
gests that the forest canopy was formed primarily by 
taxodiaceous, cupressaceous, and podocarpaceous coni- 
fers [100,119,120], a composition typical of most Late 
Cretaceous warm temperate forests [117]. Angiosperm 
shrubs may [98] or may not [121]) have formed an 
understory, alongside tree ferns and gymnosperm 
saplings [100]. The herb layer would have included ferns, 
lycopods, angiosperm herbs, and gymnosperm saplings, 
and ground cover comprised mosses, lichens, fungi, 
hornworts, and decaying vegetable matter [100]. 

Opinions vary about the degree to which habitat struc- 
ture influenced the regional distribution of the Late 
Cretaceous herbivorous dinosaurs. Some [97,101,112] 
argued that the megaherbivorous forms were likely 
restricted to feeding in open habitats, partly as a result 
of their large sizes. However, it appears that the forests 
of the Late Cretaceous were not particularly dense [96], 
and probably did not inhibit the movement of even the 
more massive herbivores [13,96]. Alternatively, Baszio 
[122] suggested that, within the DPF, ankylosaurs and 
ceratopsids occupied open habitats, whereas hadrosaurids 
lived in forested environments. His reasoning was that an- 
kylosaurs and ceratopsids, being limited in their range of 
vertical movements, could not have taken full advantage 
of stratified forest vegetation in the same way that 



Table 8 Medium-scale Mann-Whitney pairwise 
comparisons of the MAZ-2 analysis 


Ce 


Ch 


Ha (Q) 


Ha (B) 


Ce 


0.38 


0.15 


0.08 


Ch 1 




0.15 


0.08 


Ha (Q) 0.89 


0.89 




0.15 


Ha (B) 0.49 


0.49 


0.89 





N = 11, H = 8.561, p = 0.03574. 

Bonferroni corrected p-values shown in lower left triangle; uncorrected 
pairwise comparison p-values shown in upper right triangle. Significant results 
reported in bold. Abbreviations as in Tables 1 and 2. 



Table 9 Fine-scale Mann-Whitney pairwise comparisons 



of the MAZ-2 analysis 


Sty 


Cha 


Pro (Q) 


Pro (B) 


Sty 


0.77 


0.25 


0.15 


Cha 1 




0.15 


0.08 


Pro (Q) 1 


0.89 




0.15 


Pro (B) 0.89 


0.49 


0.89 





N = 10, H = 7.72,7 p = 0.052. 

Bonferroni corrected p-values shown in lower left triangle; uncorrected 
pairwise comparison p-values shown in upper right triangle. Significant results 
reported in bold. Abbreviations as in Tables 1 and 3. 
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hadrosaurids presumably could. However, it is unlikely that 
hadrosaurids could have accessed the entire forest struc- 
ture; the canopy was almost certainly out of reach, particu- 
larly if Late Cretaceous taxodiaceous and cupressaceous 
conifers grew as tall as their modern descendants (> 90 m). 
Hadrosaurids likely could forage among the shrubs of the 
forest understory, but shrubs were abundant in more open 
habitats as well [12,33,97,98,101,105]. In that case, there is 
little reason to suspect that hadrosaurids could not have 
occupied both open and closed habitats [36], alongside 
ankylosaurs and ceratopsids. 

This does not contradict the idea that certain groups 
may have preferred certain environments over others. 
Various lines of sedimentological evidence have been 
brought to bear on the matter [7,12,60,61,91,123,124]. 
There simply does not appear to have been any major 
structural obstacles to impede the movement of these 
animals. Consider that elephants, which are comparable 
in size to the megaherbivorous dinosaurs considered 
here, regularly occupy even dense forests and thickets in 
search of food [125]. In fact, their movements and feed- 
ing habits typically result in the creation of new, more 
navigable habitats [126,127]. 

Regardless of where hadrosaurids spent most of their 
time, it is likely that they usually foraged quadrupedally on 
abundant, low-lying herbage [87,89,128], occasionally 
rearing up onto their hindlimbs to feed among the angio- 
sperm shrubs. Additional evidence for bipedal feeding in 
these animals comes from the Campanian aged Blackhawk 
Formation of Utah, where hundreds of dinosaur footprints 
are preserved in association with taxodiaceous conifer and 
palm roots and fallen logs [129]. In many places, pes 
prints attributed to hadrosaurids are found straddling the 
roots. The fact that manus prints are not also found in 
these areas suggests that these animals were rearing up to 
feed on the high foliage. This bipedal feeding behaviour 
would have been particularly beneficial in instances 
where large herds of low-browsing ceratopsids were 
passing through the same area [17,20-22]. Dietary niche 
partitioning could have been achieved among hadrosaurids 
if they utilized different levels within the shrub layer, as do 
living ruminants [82]. This may also have served to limit 
niche overlap between different ontogenetic stages of the 
same species [130]. The larger feeding heights of the 
hadrosaurids suggest that these animals were able to reach 
a wider variety of plant types than other sympatric herbi- 
vores. Circumstantial evidence for diet in these animals 
comes from multiple examples of fossil gut contents 
[131-133], which preserve conifer and angiosperm browse, 
including twigs and stems, bark, seeds, leaves, and fruit. 
Probable hadrosaurid coprolites [134,135] also contain 
abundant fungally-degraded conifer wood, which would 
presumably indicate that hadrosaurids fed at ground level 
at least occasionally. However, in light of the problems 



associated with the attribution of some of these fossils 
[32,132,136,137], their interpretation as dietary residues 
must be regarded with due caution. 

Ankylosaurs, ceratopsids, and small ornithischians may 
have partitioned the herb layer by feeding height, as do 
the ungulates of the Serengeti today [78] . Ceratopsids, be- 
ing slightly taller, may have even facilitated the existence 
of the other forms by cropping the herb layer to expose 
new growth. Of course, this is a highly speculative sce- 
nario requiring further investigation. Unfortunately, no 
ceratopsid gut contents are known by which to gauge 
these ideas, but an ankylosaurid cololite from Australia is 
reported to contain fibrous tissue (probably leaves), angio- 
sperm fruits or endocarps, small seeds, and possible fern 
sporangia [138]. Ankylosaurs, it would seem, consumed 
less woody browse than hadrosaurids, which is in line with 
the interpretation given here. 

One final aspect of herbivorous dinosaur ecology bears 
consideration. Elephants are known to regularly fell trees 
up to 10 m tall to feed on the otherwise unreachable 
browse, effectively increasing the feeding envelope of these 
animals up to three times [13,125]. It is possible that the 
megaherbivores of the DPF were capable of the same be- 
haviour [13,65]. If so, tree felling may have served to in- 
crease dietary overlap between these animals, with the 
squat ankylosaurs and ceratopsids consuming foliage 
otherwise in reach of hadrosaurids alone. Unfortunately, 
while tree felling behaviour among dinosaurs is plausible, 
there is not yet any evidence supporting this speculation. 
Similarly, scenarios involving hypsilophodontids and 
pachycephalosaurids climbing trees to increase their feed- 
ing heights [122], while not impossible, are implausible 
owing to a lack of appropriate skeletal adaptation 
[139,140]. For these reasons, such highly speculative be- 
haviours are not considered further. 

Evolutionary palaeoecological implications 

There is no convincing evidence that feeding height 
stratification, as revealed by reconstructed MFH, played 
as significant a role in facilitating herbivorous dinosaur 
niche partitioning in the DPF as previously assumed 
[12,13]. Despite the 18 genera considered here— six or 
more of which typically coexisted at a time [51]— only 
four statistically distinct MFHs are detected. If niche 
partitioning did allow herbivorous dinosaurs from the 
DPF to coexist, it may have been achieved by other 
means in addition to feeding height stratification. Al- 
though this hypothesis has yet to be subjected to rigor- 
ous testing, multiple morphological features have been 
proposed to have fostered the coexistence of these herbi- 
vores. For example, Carpenter [38-41] suggested that 
differences between the tooth and beak shapes of 
ankylosaurids and nodosaurids may have allowed these 
taxa to specialize on different plant types. Similarly, 
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differences between centrosaurines and chasmosaurines in 
cranial [37], mandibular [45], and beak [8,43] morphology 
have been cited as evidence for dietary niche partitioning. 
Finally, variations in beak shape [7,34,36,44], tooth morph- 
ology [36,46], and skeletal proportions [36] are thought to 
have enabled hadrosaurines and lambeosaurines to forage 
differentially. Many of these assumptions have not been 
tested and require further examination, particularly in light 
of questions regarding the significance of intraspecific vari- 
ation and the influence of time-averaging. 

The disappearance of ankylosaurs from the upper in- 
tervals of MAZ-2 of the DPF [51] suggests the possibility 
that some change in their habitat structure caused their 
displacement. Although it is by no means obvious 
whether such a change did occur, the gradual transgres- 
sion of the Western Interior Seaway over the approxi- 
mately 1.5 Ma span of the DPF undoubtedly would have 
had some influence on the palaeoflora. It may be that some 
of the herbaceous plants preferred by the ankylosaurs 
disappeared, but this scenario is difficult to test at present. 

Overall, the distribution of herbivore MFHs changed 
minimally over the course of the DPF. Rather, MFHs 
were quite stable in spite of rapid and continual species 
turnover, and roughly the same ratio of low to high 
browsers was upheld (Figure 3B,C). This, in turn, sug- 
gests that time-averaging does not completely obscure 
palaeoecological signals within the DPF, other than to ar- 
tificially inflate estimates of standing crop biodiversity. It 
also suggests that the MFHs maintained by their respect- 
ive species were evolutionarily stable strategies, and may 
reflect correlated stability in the growth habits of the 
surrounding plants. Major changes in habitat structure 
do not appear to have occurred until the beginning of 
the Paleocene [101,105], underscoring the importance of 
low-growing herbage in sustaining Late Cretaceous 
herbivore faunas. 

Finally, it must be noted that, while differences in esti- 
mated MFHs are consistent with the hypothesis of feed- 
ing height stratification, they are not sufficient for this 
hypothesis to be true. It is possible that, despite these 
differences, all herbivore taxa from the DPF spent most 
of their time feeding at ground level [33,34]. In this 
sense, they may be compared to grazing ungulates, 
which spend most of their time feeding on low grasses, 
despite being physically capable of reaching higher 
browse. Therefore, although it seems likely from an eco- 
logical perspective that the herbivorous dinosaur fauna 
of the DPF exhibited some form of feeding height strati- 
fication, competing hypotheses about the role of this 
mechanism in the facilitation of niche partitioning are 
underdetermined [141] by the available evidence. To re- 
ject the null hypothesis of no feeding height stratifica- 
tion, it would be necessary to show that the herbivorous 
dinosaurs browsed to their full potential, utilizing their 



entire reconstructed MFHs, and did not simply spend all 
their time feeding at ground level. Unfortunately, this 
type of behaviour simply does not fossilize. Nevertheless, 
it may be possible to approximate the amount of time 
spent feeding at different heights by observing other as- 
pects of morphology. For example, it has been shown 
that low-level grazers often possess a suite of cranial 
characteristics that allow them to efficiently crop short 
grass, including wide, ventrally deflected muzzles, elong- 
ate faces, transversely wide paroccipital processes, deep 
mandibles, and tall withers [142-150]. Similarly, primates 
feeding close to the ground generally possess narrower 
dental microwear scratches than those feeding higher up in 
the forest canopy, a function of mean particle size and the 
ratio of soil particles to phytoliths [151-153]. These same 
features might be sought among the herbivorous dinosaurs 
to more accurately determine their browsing habits. 

Conclusions 

Niche partitioning between herbivorous dinosaurs from 
the DPF could have been facilitated in part by feeding 
height stratification, as is commonly observed among the 
ungulates of the Serengeti [78,82]. Small ornithischians 
(cf. Orodrorneus, Unescoceratops, pachycephalosaurids), 
ankylosaurs, and ceratopsids were generally restricted to 
feeding below 1 m, and although there is limited data 
available for the first group, ankylosaurids and ceratopsids 
may have partitioned the herb layer. Hadrosaurids may 
also have fed on herbs, but could have avoided competing 
with the other forms by feeding on taller angiosperm 
shrubs in times of duress. Sympatric hadrosaurid species, 
or ontogenetic stages, also could have avoided competing 
with one another in this way. There is, as yet, no evidence 
to suggest that any of these taxa were restricted to feeding 
in either open or closed habitats. 

Despite the evidence for feeding height stratification 
among ankylosaurs, ceratopsids, and hadrosaurids, there 
is currently no evidence for feeding height stratification 
within these clades (with the possible exception of quad- 
rupedal versus bipedal hadrosaurids described above). 
We note, however, that the problem of underdetermin- 
ation does not allow us to reject the possibility that 
some form of feeding height stratification occurred 
within the herb layer. It is also possible that the coexist- 
ence of the herbivorous dinosaurs from the DPF was 
facilitated by other means in addition to feeding height 
stratification. The prolonged stability of the palaeoeco- 
logical relationships identified here supports the conten- 
tion that the fossil assemblage from the Dinosaur Park 
Formation constitutes a chronofauna. 
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